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Abstract. We study the collective behavior of inclusions inducing local anisotropic curvatures in a flexible 
fluid membrane. The iV-body interaction energy for general anisotropic inclusions is calculated explicitly, 
including multi-body interactions. Long-range attractive interactions between inclusions are found to be 
sufficiently strong to induce aggregation. Monte Carlo simulations show a transition from compact clusters 
to aggregation on lines or circles. These results might be relevant to proteins in biological membranes or 
colloidal particles bound to surfactant membranes. 
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The interplay between structural features and iV-body 
interactions is a general physical problem arising in many 
different contexts, e.g., crystals structure 0]], magnetic 
atom clusters [|,[|, colloids in charged fluids |Q, poly- 
electrolyte condensation J5Uq] , and protein aggregation in 
biological membranes 0. TV-body interactions can some- 
times yield spectacular effects: non-pairwise summability 
of charge fluctuation forces can dramatically affect the 

| stability of polyelectrolyte bundles ||; three-body elas- 
tic interactions may induce aggregation of membrane in- 
clusions, although two-body elastic interactions are repul- 

\ sive [0. In a system able to kinetically achieve equilib- 
rium, the clusters formed are usually compact, however 

' certain interactions may favor tenuous clusters. For in- 
stance, a recent iV-body study has shown that above a 
critical strength of three-body interactions, the state of 
minimum energy is one in which all the particles are on 
a line g. It has also been observed recently that mem- 
brane mediated interactions can induce one-dimensional 
ring-like aggregates of colloidal particles bound to fluid 
vesicle membranes |J. 

Manifolds embedded in a correlated elastic medium 
can impose boundary conditions, or modify the elastic 
constants. This usually gives rise to mean-field forces, which 
are due to the elastic deformation of the medium, and to 
Casimir forces, which are due to the modification of its 
thermal fluctuations. Such interactions are generally non 
pairwise additive |1(J . The clastic interactions between de- 
fects in solids jJT[ or in liquid crystals [1^] are well known 
examples of mean-field forces. Casimir forces exist be- 
tween manifolds embedded in correlated fluids, such as 
liquid crystals and superfluids [113 
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1 OH, or critical mix- 



tures g. Another interesting eYaxnpk is the interaction 
between inclusions in flexible membranes |l6[ : it has been 
shown that cone shaped membrane inclusions experience 



both long range attractive Casimir interactions and re- 
pulsive elastic interactions falling of as R with separa- 
tion R 0. 

In this Rapid Note, following Netz fl8| , we give exact 
results concerning the long range multi-body interactions 
among membrane inclusions that break the bilayer's up- 
down symmetry. However, rather than supposing that the 
inclusions simply induce a local spontaneous curvature, 
we assume that the inclusions set a preferred curvature 
tensor |17 1!|. This model is more realistic: the "prefer- 
ence" of a conically shaped inclusion is c\ = c 2 = Co, 
where c\ and c 2 are the membrane principal curvatures, 
rather than the weaker condition c\ + c 2 = 2cq assumed 
in Rcf. |18fl . In addition, the imposed curvature tensor can 
be anisotropic, thus describing inclusions that break the 
in-plane symmetry. In a first part we calculate the exact 
Casimir and mean-field two- and three-body interactions 
between such anisotropic inclusions. Then, the collective 
behavior of identical inclusions is investigated by means of 
a Monte Carlo (MC) simulation, using the full TV-body in- 
teraction energy plus a hard-core repulsion modeling the 
simplest repulsive short-range interactions |20(| . Our re- 
sults could be relevant to understanding the aggregation 
and organization of proteins in biological membranes, or 
colloidal particles bound to surfactant membranes H . 

Let us consider a system of N anisotropic inclusions 
embedded in a flexible fluid membrane, in which they are 
free to diffuse laterally. In many situations, the surface 
tension is negligible and the membrane shape is governed 
by the Hclfrich curvature energy h n — ^n(a + c 2 ) 2 + 
R c\Ci plj , where k is the bending rigidity, and k the Gaus- 
sian modulus. For biological membranes, k ~ 30 T, while 
for surfactant membranes it can be as small as a few T (T 
will denote throughout the temperature in energy units). 
We model the membrane shape by a simple parametric 
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surface (r, u(r)), where r is a vector in the (x, y) plane and 
u(r) the normal displacement field along z. To quadratic 
order in u, the Helfrich Hamiltonian takes the form 



Ho 



r/r 



k(V 2 u) 2 + Kdet(VVu) 



Its correlation function, (u(O)u(r)) = (T/k) G(r) is given 
by the Green function 



G(r) = (V 4 )-^(r) = G(0) - ^r 2 In (?j 



(2) 



where L is a long wavelength cut-off, which is comparable 
with the size of the membrane. 

Typical membrane inclusions have a central hydropho- 
bic region spanning the hydrophobic core of the mem- 
brane, and two polar extremities protruding outside J22J. 
Assuming a strong coupling between the lipids and the 
inclusion's boundary p3| , p0[ , which in general is not cylin- 
drical, we model membrane inclusions as curvature sour- 
ces 1 17 , point-like as in p9[ . Indeed, the size of the pro- 
teins is comparable with the short wavelength cut-off, i.e., 
the membrane thickness. Even large particles can be treated 
as point-like inclusions, provided one uses the curvature 
energy coarse-grained to the size of inclusions with its cor- 
responding renormalized bending rigidity (24, 25j. Hence, 
for an inclusion located at r„, we enforce the local condi- 
tion 



VVt 



= Q„ = 



K n -J n 
K n + J n 



(3) 



where VVu is the curvature tensor of the membrane, 
formed by the second derivatives of u, and Q ra is the cur- 
vature tensor locally imposed by the inclusion. Here, Q n is 
written in the frame of reference where it is diagonal, K n 
being the mean-curvature and J n the anisotropic curva- 
ture of the inclusion. While isotropic inclusions correspond 
to J n — 0, "hyperbolic" and wedge-shaped inclusions cor- 
respond to K n = and K n = J„, respectively. In the fol- 
lowing, we may assume J n > without loss of generality. 
We define the orientation of inclusion n by the direction 
of the lowest principal curvature, i.e., K n — J n . The latter 
is defined modulo a rotation of angle n. Indeed, even if 
the inclusion bears an in-plane polarity, its coupling with 
the membrane curvature is apolar by symmetry. 

Within this model, we can calculate the free energy T 
of a membrane with N inclusions in a non perturbative 
way. Introducing an external field h(r), eventually set to 
zero, the free energy is given by 



exp[-— ] = / £>[w]exp 



±Jdvu{nV 4 



(4) 



®ln 




Fig. 1. Definition of the mutual orientations of three inclu- 
sions labeled n,m, and I. The angle 9 n i is the orientation of 
inclusion n with respect to the line joining inclusions n and I, 
etc. All angles are defined modulo it, since any orientation and 
its reverse are equivalent. 



This integral will be regularized by introducing a micro- 
scopic cutoff a of the size of the membrane thickness. The 
Gaussian curvature energy, which yields a constant contri- 
bution depending on the boundary conditions at the edge 
of the membrane and a constant self-energy per inclusion 
~ -Tra 2 R(K% - J 2 ) (consistent with Ref. Q), has been 
discarded. 

Replacing the delta functions by their Fourier repre- 
sentation (for details on this method, see, e.g., Ref. |10[), 
and integrating out the field w(r), we obtain 



T[h] = r ^ln{dctM) + ^Kb'M- 1 b' t , 



b' = b+ fdrh(r)c{r). 



(G) 

(7) 



M is a 37V x 3 AT matrix formed by the blocks M„ p = 
D'DG(r n - r p ), where D = (dl,d x d y ,dy); b and c(r) 
are 3N vectors with blocks b n = [(Q„)n, (Q„)i2, (Qn)22] 
and c„ = D G(r — r„), respectively. 

The average shape of the membrane is given by 



Mr)) 



6F[h] 



Sh(r) 



cfrJM^b' 



(8) 



h=0 



Although not in explicit form, Eq. (|6|) allows to calculate 
exactly the W-body interaction, by inverting the matrix 
M whose elements are simple functions of the distances 
between particles. M and b do not depend on T, therefore 
the first term in (|6|) corresponds to the fluctuation induced 
Casimir interaction, while the second term represents the 
mean-field elastic interaction. 

Expanding the Casimir interaction T G = ^Tln(det M) 
to fourth order in l/r np , where r np = |r„ — r p |, yields 



in which the constraints are implemented by delta func- 
tions in the measure 



iV 



V[u] - V[u] [] 6 (VVu| rn - Q n ) . 



(5) 



T? 



-3TV — , 



(9) 



n , p 



in which the prime indicates restriction of the sum to 
different inclusions. This fourth-order Casimir interaction 
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is pairwise additive and corresponds exactly to the sum 
of the two-body interactions obtained in Refs. [0@ f° r 
isotropic inclusions. It does not depend on the locally im- 
posed curvature nor on the orientation of the inclusions 
(unlike the Casimir interaction between two rod-like inclu- 
sions pq| ) . This interaction is quite weak, since for r > 2a 
it is only a fraction of T and thus it is overwhelmed by 
diffusion. 

Setting h — in the second term of (g) yields the 
mean- field elastic energy T mt = ^KbM" 1 ^, which arises 
from the bending deformation of the membrane. To second 
order in l/r np , we obtain (see Fig. [I] for the definition of 
the angles) 



•rmf 

• / "2 — 



,' a 4 

\ZJ n Jp COSyZUpfi ZUjipj 



+K n J p cos 29 pn + K p J n cos 29 np 



(10) 



At this order, the mean-field interaction is pairwise ad- 
ditive. It depends on the imposed curvature tensors and 
on the orientations of the principal curvatures. To fourth 
order, the mean-field interaction includes both two-body 
interactions: 



T 
j t 



mi' 

42 



' a" 



cos49 np )J* 



(5 + cos46» p „)^ 2 



n,p n P 

K 2 n + Kl+AK n J n cos29 np + 4K p J p cos20 pn ] , (11) 



and three-body interactions: 



? 43 = -T K Yl ^^-UnJp[4:COs(29r l 

6 n, P ,e r ^ r tp 



-26„p+2a 



rip) 



2 cos 26 n i cos 29 p e\ + 2K n J p coa(29 p i+2a^ p ) 
+ 2K p J n cos(26 ni +2ai p )+K n K p cos2ai p \ + peim(12) 

These terms can be attractive or repulsive, depending on 
the curvatures and orientations of the inclusions. Setting 
J n — J P = in ( pd| ) , we recover the result of Refs. [ p"7||l9[ ] for 
two isotropic inclusions; however, we also notice that the 
three-body interactions arc in general of the same order 
as the two-body interactions. 

Let us study the lowest-order mean-field interaction 
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between two identical anisotropic inclusions. Setting 
= K-2 = K and J\ = J2 = J, we may assume K > and 
J > without loss of generality. UK ^ 0, the energy is 
minimal when #12 = #21 = 0. The inclusions therefore tend 
to align their axis of smallest principal curvature (smallest 
in modulus) parallel to their separation vector. If K — 
0, the energy is minimal whenever 9\2 = 9%\, In both 
cases the interaction is attractive, contrary to the repulsive 
interaction found between two isotropic inclusions |17|Jl^|. 
If #12 = and 621 = is arbitrary, the interaction is 
repulsive for 6* < \6\ < it -9*, where cos 29* = —K/(K + 
2 J), and otherwise attractive; the relation 7r/4 < 9* < 
tt/2 holds, where the lower limit corresponds to K — 0, 
and the upper limit to J — > 0. When #12 = 621 = 9, 
the interaction is also repulsive for 9 C < \9\ < it — 9 C 



a) 



b) 



c) 



d) 




Fig. 2. Typical equilibrium aggregates obtained from MC sim- 
ulation of iV = 20 identical inclusions. The corresponding val- 
ues of (K, J) are indicated in the phase diagram of Fig. 0. The 
mesh size corresponds to the microscopic cut-off a. The bars 
indicate the orientations of the inclusions. 



and attractive otherwise, with cos26> c = —J/K. Again, we 
have 7r/4 < 9 C < tt/2, where the lower limit corresponds 
to J — > and the upper limit to J = K (wedge-shaped 
inclusions). If J > K the interaction is always attractive 
when 9\2 = #21 • 

As seen above, the lowest-order interaction between 
two identical anisotropic inclusions is already quite com- 
plex, however it is essentially attractive. Hence, one ex- 
pects aggregation of the particles when the interaction 
energy overcomes the entropy of mixing. We noticed that 
the case K — is special, since the minimum energy is de- 
generate. Thus, one might expect compact clusters in this 
case, but for larger values of K it should become increas- 
ingly favorable to orient the inclusions along a common 
line joining them. To check these ideas, we have performed 
MC simulations on finite samples of N = 20 identical in- 
clusions. The positions and directions of principal curva- 
tures were varied continuously, in order to avoid lattice 
artifacts. The energy was calculated from the iV-body in- 
teraction (o), with h = 0, by numerically inverting the 
3N x 3iV matrix M. When inclusions are close to contact, 
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Fig. 3. Phase diagram of identical anisotropic inclusions ob- 
tained by MC simulations for different values of the inclusions 
mean-curvature K and anisotropic curvature J. (G) gas phase, 
(A) compact aggregates, and (P) linear polymer-like aggre- 
gates. In the dashed region, we found branched and ring-like 
aggregates, or small clusters close to the gas phase. 



additional short-r ang e interactions should be taken into 
account [|l6|j2C| , |23| , P7| . To mimic them, we have added a 
repulsive hard-core potential (diameter 4a), in agreement 
with the results of Ref. |2(J in which a short-range energy 
barrier was obtained from a microscopic model. Obviously, 
in the case of attractive short-range interactions, we can- 
not predict the aggregate shapes within the present model. 
The results of our MC studies merely give the aggregation 
tendencies resulting from the long-range interactions. 

We started the simulations with a low-density distri- 
bution of inclusions having both random positions and ori- 
entations, then we let the system equilibrate using ~ 10 7 
MC steps. The relevant dimensionless parameters are K = 
(n/TY^Ka and J = (k/T) 1 / 2 Jo. For weak anisotropy, 
i.e., J < 1, the interaction was either repulsive or too 
weakly attractive to induce aggregation; for J > 1, the 
interaction was sufficiently attractive for the inclusions to 
aggregate. For small K, we obtained compact clusters of 
inclusions, inducing some kind of decorated "egg-carton" 
structure of the membrane (Fig. ||a). As K increased, we 
found a transition from compact aggregates to polymer- 
like ones (Fig. 0d). Remarkably, we also found ring-like 
aggregates, featuring a budding mechanism (Fig. gb,c). 
The boundaries between the different aggregates shapes 
are depicted in the phase diagram of Fig. 0. Performing 
simulations with N — 60, wc observed relatively weak fluc- 
tuations of the linear aggregates, indicating a persistence 
length of several hundred of inclusions or more. This sug- 
gest that "polymers" made of wedge-shaped inclusions can 
acquire a persistence length in the micron range. Long- 
range order might exist in the form of parallel inclusion 
lines, since for a two dimensional object, the 1/r 2 inter- 
action is marginally long-range. In all these aggregates, 
the multi-body interactions are of the same order as the 
pairwise interactions, but the Casimir interaction is not 
determinant for the types of aggregates obtained. 



Our MC simulations concerned only systems of identi- 
cal inclusions, however they may easily be extended to the 
case of inclusions inducing different curvatures. It would 
also be interesting to study the effect of the long-range 
anisotropic interactions in the situation of kinetically con- 
strained aggregation. In the case of a film with surface 
tension, we found similar anisotropic interactions, but of 
shorter range ~l/r 4 at lowest order pqj . 

We thank R. Bruinsma and J. Prost for useful discus- 
sions. P. G. D. was supported by the Research Council of 
Norway. 
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